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1 Introduction 


In Part I of this series of papers [TOlimfT^ (partly inspired by [9]) we have considered 
the general framework of local (and global) zeta regnlarization of a nentral scalar 
held in a d-dimensional spatial domain fl, in the environment of {d + l)-dimensional 
Minkowski spacetime; the possible presence of an external potential V was taken 
into acconnt as well. 

In the present Part III we consider a massless held on = R'^, and choose for V a 
harmonic potential. The potential is assumed to be isotropic, i.e., to be proportional 
to the squared radius |xp (x E R'^); nevertheless, our approach could be extended 
with little ehort to anisotropic harmonic potentials and to the case of a massive 
scalar held. 

The main result of this paper is the renormalization of the vacuum expectation value 
(VEV) for the stress-energy tensor, obtained applying the general framework of Part 
I to the present conhguration; we also discuss the total energy in this conhguration. 
After producing general integral representations for both the renormalized stress- 
energy VEV and the total energy, we consider in detail the cases of spatial dimension 
de {1,2,3}. 

The above mentioned integral representations are derived analytically and are fully 
explicit; however, the integrals therein must be computed numerically. We have 
performed these latter computations for d E {1,2,3}, using Mathematica. 

Using the previously cited integral representation, it is also possible to derive the 
asymptotics for the stress-energy tensor components when the radius |x| goes to 
zero or to inhnity; we derive, as well, remainder estimates for these asymptotic 
expansions as well. 

The idea to replace sharp boundaries with suitable background potentials is well- 
known in the literature on the Casimir effect. Typically (see, e.g., PEI [151 Ellin]), 
delta-like potentials are introduced in order to mimic boundary conditions in a 
“physically more realistic” framework; the ultimate purpose is to obtain less singu¬ 
lar behaviours of the renormalized quantities, avoiding, e.g., boundary divergences 
such as the ones pointed out in Parts I and II [TUlfTTj . The case of a scalar held inter¬ 
acting with an external harmonic potential has been formerly considered by Actor 
and Bender BE], who have determined the renormalized VEV of the total en¬ 
ergy via global zeta regularization; to the best of our knowledge, local aspects such 
as the stress-energy tensor have never been considered previously for the present 
conhguration. 

The paper is organized as follows. In Section [2] we present a summary of basic results 
from Part I; as in Part II, the purpose of the summary is to save the reader from 
the tedious task of recovering from Part I the general identities applied here. In 
Section |3] (and in the related Appendix [^ we discuss the general approach to treat 
the case of a scalar held interacting with a classical, isotropic harmonic potential in 
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arbitrary spatial dimension d] this includes general expressions for the renormalized 
stress-energy VEV, its asymptotics for |x| small or large and for the renormalized 
VEV of the total energy (more precisely, for the bulk energy as dehned in Part I). 
The main ingredients in the derivation of the above results are: 

i) the general results of Part I relating the zeta regularization to the heat kernel of 
fundamental operator A = — A -|- E(x); 

ii) the fact that, when E(x) is harmonic, the related heat kernel is the well-known 
Mehler kernel [18] (also see [H iGl 181 fT 6 ]h 

In Section 01 the analysis presented in the preceeding sections is specihed to the 
cases of spatial dimension d G {1, 2, 3}. 

We point out that results on the renormalized bulk energy agree with the ones of 
Actor and Bender |2| for d G {1,2,3}. Nevertheless, it should be mentioned that 
these authors compute the analytic continuations required by the zeta approach 
using ad hoc results on continuation the special functions involved in this specihc 
case; on the contrary, here we are applying mechanically the general setting discussed 
in Part I, in the spirit of the present series of papers. 


2 A summary of results from Part I 

2.1 General setting. Throughout the paper we use natural units, so that 

c=l, h=l. (2.1) 

Our approach works in (d-l-l)-dimensional Minkowski spacetime, which is identihed 
with using a set of inertial coordinates 

X = (x'")^=o,i,...,d = (x°,x) = (t,x) ; (2.2) 

the Minkowski metric is {rjnu) = diag(—1,1,..., 1). We fix a spatial domain O C R'^ 
and a background static potential E: O — R. We consider a quantized neutral, scalar 

held 0 : R X O —)■ Asa(i?) (i? is the Fock space and Csa{d) are the selfadjoint operators 

on it); suitable boundary conditions are prescribed on dfl. The held equation reads 

0 = (-9ttA - l/(x))0(x,f) (2.3) 

(A := i® i'h® d-dimensional Laplacian). We put 

A:=-A + V, (2.4) 

keeping into account the boundary conditions on dfl, and consider the Hilbert space 
L^(f 2 ) with inner product {f\g) ■= d'x f {-x)g{-x). We assume A to be selfadjoint 
in L^(f2) and strictly positive (i.e., with spectrum a{A) C [e^, +oo) for some £ > 0). 
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We often refer to a complete orthonormal set {Fk)k£ic of (proper or improper) eigen- 
fnnctions of A with eigenvalnes {ul)k£ic for all kEJC). Thus 

Ffc : — )■ C; AFk = uj^Fk ; ^ 

{Fk\Fh) = S{k, h) for all fc, h G /C . 

The labels fc G /C can include both discrete and continuous parameters; dk indi¬ 
cates summation over all labels and 6{h, k) is the Dirac delta function on /C. 

We expand the held 0 in terms of destruction and creation operators corresponding 
to the above eigenfunctions, and assume the canonical commutation relations; |0) G 
5 is the vacuum state and VEV stands for “vacuum expectation value”. 

The quantized stress-energy tensor reads G R is a parameter) 

:= (1 - 2^)d^$od^$- Q - 2^^ - 2^ 0 o ; (2.6) 

in the above we put A o B := {1/2){AB -|- BA) for all A,B G Fgaid), and all the 
bilinear terms in the held are evaluated on the diagonal (e.g., d^cj) o d^ct) indicates 
the map x i-G- d^(j){x) o du(l>{x)). The VEV (0|T^i,|0) is typically divergent. 

2.2 Zeta regularization. The zeta-regularized Reid operator is 

0“ ;= (K-2>f)-V40 ^ (2.7) 

where A is the operator fl2.4|] . u G C and k > 0 is a “mass scale” parameter; note 
that 0“|ti=o = 0, at least formally. The zeta regularized stress-energy tensor is 

A - {i-20s,roay'‘-A-2(),,^,(ay'‘axr+v{rf)-2(rod,,d'‘. ( 2 . 8 ) 

The VEV (0|T“^|0) is well dehned for lHu large enough (see the forthcoming subsec¬ 
tion |23]); moreover, in the region of dehnition it is an analytic function of u. The 
same can be said of many related observables (including global objects, such as the 
total energy VEV). 

For any one of these observables, let us denote with F{u) its zeta-regularized version 
and assume this to be analytic for u in a suitable domain Uq G C. The zeta approach 
to renormalization can be formulated in two versions. 

i) Restricted version. Assume the map Wq —)■ C, u F(u) to admit an analytic 
continuation (indicated with the same notation) to an open subset U C C with 
7/^0; then we dehne the renormalized observables as 

Fren := -^(0) • (2.9) 
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ii) Extended version. Assume that there exists an open subset UdC with Uq <zU, 
such that OgW and the map u^Uq i—)■ iF{u) has an analytic continuation to W\{0} 
(still denoted with iF). Starting from the Laurent expansion J^{u) = Ylik=-oo^kU^i 
we introduce the regular part {RP J^){u) := dehne 

Rren-.= {RPJ^m . ( 2 . 10 ) 


Of course, if T is regular at u = 0 the defnitions (I2.9p (I2.inp coincide. 

In the case of the stress-energy VEV, the prescriptions (i) and (ii) read, respectively. 


(0|f„„(a:)|0)„„ := (0|f"„(x)|0) , 


( 2 . 11 ) 


{0\%..(x)\0)ren BP 


/0|f“ (l)|0) , 


( 2 . 12 ) 


2.3 Conformal and non-conformal parts of the stress-energy VEV. These 
are indicated by the superscripts (0) and (■), respectively; they are dehned by 


m^u\^)ren = + (^-^d) (010),en , 


where we are considering for the parameter ^ the critical value 


id-.= 


d — l 


(2.13) 


(2.14) 


2.4 Integral kernels. If S is a linear operator in L^(0), its integral kernel is the 
(generalized) function (x,y) e OxO i—)■ B{x,y) := (5x|i3(5y) ((5x is the Dirac delta 
at x). The trace of B, assuming it exists, fulhlls TrB = f^dxB(x,x). 

In the following subsections .4, is a strictly positive selfadjoint operator in L^(Q), with 
a complete orthonormal set of eigenfunctions as in Eq. fl2.5p . In typical applications, 
A is the operator fl2.4p . 

2.5 The Dirichlet kernel and its relations with the stress-energy VEV. 

For (suitable) s G C, the s-th Dirichlet kernel of A is 

r dh _ 

D,(x,y) := ^■"(x,y) = / — Fk{x)Fk{y) . (2.15) 

Jk 

A = —A -|- V (with V a smooth potential) is strictly positive, the map Fs( , )■ 
D X D —)■ C, (x, y) i-G- Ds(x, y) is continuous along with all its partial derivatives up 
to order j G N, for all s G C with 3fJs > d/2 + j/2. 

Recalling Eq. fl2.8p . the regularized stress-energy VEV can be expressed as follows: 


{0\T^,{x)\0) = n^ 


+? C!i=i(x,y)+ --q(5“S/+r(x))r>s«(x,y) 


(2.16) 


y=x 
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= K 


u 



(0|To“(x)|0) = (0|7;“(x)|0)=0, (2.17) 

(0|7;“(x)|0) = (0|f™(x)|0) = 

-i)5ij{p^{^,y) - {d^'dy, + V{^))D^{^,y)^ + ^2 13 ) 


((0|T“j,(x)|0) is short for (0|T“^(t,x)|0); indeed, the VEV does not depend on t); 
of course, the map —)■ C, x 1 —)■ (0|T“^(x)|0) possesses the same regularity as 
the functions x G ff e-)■ (x, x), x) {z^w any two spatial variables); 

so, due to the previously mentioned results, x 1 —)■ (0|T“^,(x)|0) is continuous for 

> d + 1. 

The above framework relates the renormalized stress-energy VEV (0|T^i,(x)|0)ren : = 
7?P|u=o{0|T“j,(x)|0) to the analytic continuation at n = 0 of the maps u 1 —)■ (x, y) 

and of their derivatives. 

In the sequel we will also consider the total energy VEV and express it in terms of 
the trace Tr^“®, fulhlling 


Tr^ * = / (ixZls(x,x) . (2-19) 

Jn 

2.6 The heat kernel. For t G [0, + 00 ), this is given by 

77(t;x,y) := e"‘^(x,y) = [ d/c Ffc(x)Ffc(y) . (2.20) 

Jic 

If .4, = —A -|- V (V smooth) is strictly positive, the map d^(t; , ) : fl x fl —)■ C, 
(x, y) I— )■ A(t;x, y) is continuous along with all its partial derivatives of any order, 
for all t > 0 . 


The heat trace (assuming it to exists) is 

h:(t) := Tre“^ = / dxh:(t;x,x). (2.21) 

Jn 

2.7 The Dirichlet kernel as Mellin transform of the heat kernel. For suit¬ 
able values of s G C (see Part I), there holds 

1 ^+00 

^^(x,y) = :p^/ dtP"^h:(t;x,y) . (2.22) 

r(s) Jo 

Similar results hold for Tr 4.“^; for example, using the heat trace K{t) of Eq. (I2.2ip . 
we obtain 

1 p+00 

TtA-^ = -— dtp-iiP(t). (2.23) 

r(s) Jo 
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2.8 Analytic continuation of Mellin transforms via integration by parts. 

Let J- ; (0, +oo) —)■ C be a function of the form 

^(t) = iH(t) (2.24) 

for some p G C and some smooth function "H : [0, +cx3) —)■ C, vanishing exponentially 
for t —?• +CX3. Consider the Mellin transform of i.e. 


fm(cT) 


r *+00 


di J'(t) ; 


(2.25) 


this is an analytic function of a for JJcr > 3?p. Integrating the above equation by 
parts n times (for any n G {1, 2, 3,...}) we obtain the relation 


9Jl(a) = 


(- 1 )’ 


r"+oo 


{a-p)...{a-p+n-l) Jq 


rjnqj 

(2.26) 


giving the analytic continuation of DJl{a) to a meromorphic function of a in the 
region — n}, possibly with simple poles at a G {p, p — 1, , p — n + 1}. 

The above results can be emplc^ed to obtain the analytic continuations of the reg¬ 
ularized stress-energy VEV (0|T“j,(x)|0) (treating x G as a hxed parameter) and 
of the trace Tr for example, assuming the heat trace K{t) to have the form 

/nt) = -t H{i) , (2.27) 

for some p G R and some smooth function H : [0, -|-oo) — )■ R, vanishing rapidly for 
t —)■ -|-oo, starting from Eq. fl2.23p we obtain (for n G {1, 2, 3,...}) 


TrA"" = 


-ir 


p+oo 


T{s){s-p)...{s-p+n-l) Jq 




rjn TJ 
-p+n-1 ^ ^ 

di^ 


(t). 


(2.28) 


2.9 The case of product domains. Factorization of the heat kernel. Let 

A'.= —A + V and consider the case where 

X 3 X = (xi, xa), y = (yi, ya) , (2.29) 

l/(x) = l/i(xi) + l/ 2 (xa) (2.30) 

(fia C R‘^“ is an open subset, for a G {1,2}; di + d 2 = d); assume the boundary 
conditions on dVL to arise from suitable boundary conditions prescribed separately 
on dfli and dfla so that, for a = 1, 2, the operators 

Aa :=-A, + E(x„) (2.31) 
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(with Aa the Laplacian on (]„) are selfadjoint and strictly positive in L‘^{Qa)- Then, 
the Hilbert space L‘^{Q) and the operator A can be represented as 

(g) , A = Ai(S)1 + 1(^A2. (2.32) 

This implies, amongst else, that the heat kernels i^(t; x, y) := e^(x, y), ii'(j(t; x^, y^) 
:= e‘-^“(xa,ya) (a = 1,2) are related by 

K(t;x,y) = it'i(t;xi,yi)it' 2 (t;xi,X 2 ) . (2.33) 

Similarly, writing K{t), Ka{t) (a G {1,2}) for the heat traces of A and Aa {a G 
{1,2}), respectively, we have 


K{t)=K,{t)K2{t) . (2.34) 

Let ns briefly recall that a special snbcase of the present framework is the case of a 
slab, i.e., 

H = Hi X R‘^2 3 X = (xi, xs), y = (yi, y 2 ) , V{x.) = V (xi) (2.35) 

with Hi C an open snbset {di+d 2 = d); see Part I for more details on this topic. 
2.10 The total energy. The zeta-regularized total energy is 

^“:= /'dx(0|fo“(x)|0) + ; (2.36) 

Jn 

the second eqnality is proved after dehning the regularized bulk and boundary en¬ 
ergies, which are 


:= 




K 


dx (x, x) = — Tr .4. 

2 ^ ^ 2 


\ — u 
2 


:= (i - e) £<ia(x) ^fl=±i(x.y) 


y=x 


(2.37) 

(2.38) 


If H is nnbonnded, the integral over dkl in Eq. fl2.38p has to be intended as 
lim£_,.+oo da, where (H£)£=o,i, 2 ,... is a seqnence of bonnded snbdomains with 
H^ C H^+i for i G {0,1, 2,..}, and U^H^ = H. 

Assnming E'^ and (see Eq.s fl2.37p fl2.38p l to be finite and analytic for snitable 
M G C, we define the renormalized energies by the generalized (or restricted) zeta 
approach; for example, we pnt 


Eu 

u=0 


(^or E'^^^ := 

8 


R=0 


Even j^p 


(2.39) 












2.11 Curvilinear coordinates. In some applications it is natural to employ on 
12 some set of curvilinear coordinates = q, inducing a set of spacetime 

coordinates q = (g^) = (f, q); the spatial and space-time line elements are, respec¬ 
tively, 


= aij{q)dq^dq^ ; = —dt'^ + di'^ = dq^dq'" , 

doo := -1 , dio = doi := 0 , gij{q) := ay(q) for ij G {1, ...,d} . 

The analogue of Eq. fl2.8p in the coordinate system (g^) is 


(2.40) 


T“ := (l-2Od^0“od,0“- + -2e0“oV^.0“ (2.41) 


(V^ is the covariant derivative induced by the metric fl2.40p h For any scalar function 
/ there hold (y^^- are the Christoffel symbols for the spatial metric (ajj(q)), Di the 
corresponding covariant derivative) 


V^/ = d,f , V.,/ = A,/ = d,,f - 7^4/ , 

Voi/ = doidj) = d,{dof) = V,o/ , Voo/ = doof • 


(2.42) 


In the present paper we often work in a curvilinear coordinate system; more pre¬ 
cisely we use a set of (rescaled) spherical coordinates, htting the symmetries of the 
conhguration under analysis. Most of the results of the previous subsections are 
readily rephrased in this framework. 


3 The case of a harmonic potential 

3.1 Introducing the problem. We consider the case of a massless scalar held 
on in presence of a classical isotropic harmonic potential. More precisely, we 
assume 

f2:=R^ V{x.):=k^\^W (3.1) 

where k > 0 and |x| := 'y/(adA'TTTA'(a?p ; the constant k is, dimensionally, a mass 
(or an inverse length) like the parameter k employed for the held regularization fl2.7p . 
All the considerations reported in the following could be generalized, with some 
calculational ehort, to the anisotropic case 14(x) = where ki ^ 0 for 

all i G {l,...,d}, also including slab cases where some of the ki vanish; a further 
variation would concern a massive held in one of the above mentioned conhgurations, 
so that, e.g., in place of fl3.ip we would have 14(x) := -|- /c^|xp . None of these 

generalizations will be considered, for the sake of computational simplicity. 
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3.2 The heat kernel. Even though the conhguration fl3.ip is patently spherically 
symmetric, we choose to postpone the use of spherical coordinates to the next sub¬ 
section; here we work in standard Cartesian coordinates for reasons that 

will soon become apparent. 

First of all, notice that fl3.ip is a conhguration of the product type considered in 
subsection 12.91 in particular, with a trivial generalization of the results therein, we 
have for the Hilbert space of the system and for the fundamental operator A := 
—A H, respectively, the following representations (compare with Eq. fl2.32p h 


L^(R'^) = (g)L^(R) , 

i=l 


( 3 . 2 ) 




1(8.4,i , .4,1 := —^ + on h^(R) . (3.3) 

dx^ 


A = Ai®l®...®l + + 

In this situation, the heat kernel of A is given by (compare with Eq. fl2.33p i 

d 




( 3 . 4 ) 


2=1 


where Ki is the heat kernel of Ai, the latter is the familiar Mehler kernel, that is 


Ki{t;x,y) = 


k 


Ay27rsinh(2fc2t) 


exp 




x‘^+ y'^ 


X 


y\ 


2tanh(2fc^t) sinh(2fc2t) / 


( 3 . 5 ) 


(see [0], using the equations therein with the replacement (t,x,y) —>■ {kH, kx, ky)). 
The last two relations imply (with x ■ y := 


^(t;x,y) = 


k 


sinh(2/c2t) 


exp 


|xr+ |y 


x-y 


tanh(2fc2t) sinh(2/c2t) / 


( 3 . 6 ) 


A similar analysis can be performed for the heat trace. Firstly note that (compare 
with Eq. fl2.34p i 

d 




2=1 


( 3 . 7 ) 


where A'i(t) indicates the trace for the reduced one-dimensional problem; for the 
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latter quantity we can easily deduce the expression 


^i(t) 


1 

2sinh(fc2t) 


(3.8) 


Summing up, for the heat trace of the complete d-dimensional problem we have 


K{t) = 


( - - - 

\2sinh(/c2t) 



(3.9) 


In the following subsections we explain how to construct the analytic continuation 
of the zeta-regularized VEV of the stress-energy tensor, using the expression fl3.6p 
for the heat kernel and the integral representation 02 . 221 ) for the Dirichlet kernel. 
Moreover, we show how to determine asymptotic expressions for the renormalized 
stress-energy VEV in the limit of small and large |x|, respectively. Moving along the 
same lines, we also derive the renormalized VEV of the total energy; to this purpose, 
we consider the representation 02.23|) for the trace TrM.“^ and the expression 03.9p 
for the heat trace. 

3.3 Spherical coordinates. As anticipated in subsection 13.21 we now pass to 
a set of curvilinear coordinates which best fit the symmetries of the problem (see 
subsection 12 . 111 ) : namely, we introduce the spherical “/c-rescaled” coordinates 


X I— )■ q(x) 


(r(x), 6 'i(x),..., 6 'd_ 2 (x), 9d-i 


(x)) G (0, -1 -cxd) X (0, tt) X ... X (0, vr) x (0, 27r) 

(3.10) 


whose inverse map q i—)■ x(q) is described by the equations 

kx^ = rcos{6i) , 
kx"^ = r sin( 6 'i) cos( 6 ' 2 ) , 


kx'^ ^ = rsm{6i)...sm{6d-2) cos{6d-i) , 

= r sin( 6 *i)... sin( 6 *(i_i) . 


(3.11) 


Needless to say, for d = 1 (when the system is invariant under the parity symmetry 
x^ —x^) we set 

r := k\x^\ . (3.12) 

^In order to derive Eq. (13.81) we can proceed, for example, as explained hereafter. The second 
equality in Eq. (12.211) and the explicit expression (13.51) for the one-dimensional heat kernel (along 
with some trivial trigonometric identities) allow us to infer 


^i(t) 


y^27r sinh(2fc2t) 


tanh(fe^t) 


'R 


then, Eq. (13.81) follows evaluating the above elementary Gaussian integral. 
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Let us stress that, for any spatial dimension d, we have 


r = /c lx 


(3.13) 


thus, the coordinate r is an adimensional radius. In order to avoid cumbersome 
notations, given a function hi —)■ y, x i—)■ /(x) (with Y any set), we indicate the 
composition q i—)■ /(x(q)) as q i—)■ /(q); we will use similar conventions for functions 
on n X n (see, e.g., what follows). 

Let us now consider the Dirichlet function Ds{'x,y) and the heat kernel iL(t;x,y); 


let 



q = (r,9i,. 

.., 9d-i) = (r,d) , p = (r', 9[, ..., d^.^) = (r, 0') 

(3.14) 

and put 

T ■= k‘^t e ( 0 , -fcx)) . 

(3.15) 


In the sequel we write Ds(q, p) and iL(r ;q,p), respectively, for the Dirichlet and 
heat kernels at two points x, y of (rescaled) spherical coordinates q, p, and with r 
related to t by Eq. fld.lSp . Eq. fl3.6p implies 




-^27rsinh(2r) 



/ r^+r'^ _rr' S{e)S{e')Y 

y2tanh(2r) sinh(2r) / 


(3.16) 


where S{6) and S{6') are the products of cosines and sines of the angular coordinates 
( 6 * 1 ,..., 6 *di) and {9[,9'^J of Eq. (I3.14p . corresponding to the scalar product x ■ y. 
From Eq. (12.221) it follows that 

^+oo 

Ds{q,p) = ^^ dr r"“^iL(r;q,p) ; (3.17) 

r(s) Jo 

substituting this relation into the analogues of Eq.s 02.16112.1^ for curvilinear coor¬ 
dinates (see subsection 12 .lip , we get 


^ Ld+l r+°° 

(0|T;^(q)|0)= (-j dr 77(“Hr ; q) (/i, i/G {0, r, di,..., drf_i}) (3.18) 

where the coefficients 'Hj)t^(r;q) are as follows (for i,j,hY ^ 6 *i,..., here 

D is the spatial covariant derivative) 


(^;q) := 


r 




(^Y\q)Dqhpe + Y'J 


^(^;q,p), 

p=q 


(3.19) 


Wo“’(r;q) = W';’(r;q) — 0 


(u) 


(3.20) 
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r 


+ 


(3.21) 






'wSjV;q) = 'H5?(^;q) := 


i\:(r ;q,p) . 

p=q 


Here and in the reaminder of the paper, we are implicitly nnderstanding the de¬ 
pendence on the parameter ^ for simplicity of notation: so, (r; q) stands for 

^&V;q;0- 

In the following snbsection we point ont some notable properties of the coefficients 
(r ; q); these properties will be employed in the following snbsection 13.51 in order 
to evalnate the renormalized stress-energy VEV. 


3.4 Properties of the coefficients Consider the explicit expressions (13.191 

I3.2ip for these coefficients. Using as well Eq. fl3.16p for the heat kernel expressed in 
rescaled, spherical coordinates, one can prove the following facts. 

i) For fixed r, q and any /x, i/ G {0, r, ^i,..., Od-i}-, the map u i—)■ (r ; q) is affine. 

ii) For fixed q and any n G C, /x, ix G {0, r, 6*i,..., 6'd-i}, the map r i—)■ 'Hp’t^(r;q) is 
smooth on [0, -|-oo) and exponentially vanishing for r —)■ -|-cxo. 

iii) The final expressions for the coefficients do not depend on the parameter k 
(even thongh it appears in the right-hand sides of Eq.s (13.19113.2ip h 

iv) There hold 


=0 for /i 7^ 1/ , 

^ (3.22) 

^£ 20^-2 = - = sin2(^-i-2)...sin2(ffi)-H(“^^ . 

v) For /i = IX G {0, r, ^i}, there hold 

(t i q) = ““ W (t; r) (3.23) 

where Ai‘inl{T-.r) is a polynomial In r.u of degree 1 in both these variables, with 
coefficients depending smoothly on r. 

All the above statements i)-v) conld be proved for arbitrary d, bnt we will limit 
ourself to check them in the cases d G {1,2,3} considered hereafter, for which the 
explicit expressions of all coefficients 'H|i“^(r;q) (/x = ix G {0,r, di}) will be given: 
see the subsequent subsections 14.1114.21 and 14.31 respectively. 

Before moving on, let us make a few remarks. First, consider the integral represen¬ 
tation (I3.18P for the regularized stress-energy VEV; in consequence of point iii), the 
latter VEV only depends on the parameter k through the multiplicative coefficient 
kd+l 

(/t//c)“ in front of the integral in the cited equation. Next, we note that Eq. 
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fl3.22p of point iv) indicates that the VEV (0|T“j^(q)|0) is diagonal and that the only 
independent components are those with /i = i/ G {0, r, 6i}] hnally, these components 
only depend on the radial coordinate r (and not on the angular ones {9i,9d-i}), 
in accordance with the spherical symmetry of the conhguration under analysis. 

3.5 Analytic continuation. Let us move on to determine the analytic continu¬ 
ation of the regularized VEV fid.isp . To this purpose, we refer to the framework of 
subsection I2.8| using Eq.s (12.24112.2^ with "H = ( ;x), p = {d + 3 — u)/2, we 

obtain 


pd+i 

r(iS) 


■ 1 )’ 


' u—d—3 




u—d—S 


{0|T“ (q)|0) = 

r"-|-oo 

dr t' 


-n 


k) I 


-d-3 




(3.24) 


Notice that, due to the features of (r; q) pointed out in the previous subsection, 
the integral in the above expression converges for 


> d -|- 1 — 2n , 


(3.25) 


so that Eq. fl3.24l) yields the required analytic continuation of (0|T“j,(q)|0) to the 
very same region. Since we are interested in evaluating (the regular part of) this 
analytic continuation in u = 0, we choose n so that Eq. (I3.25p holds for u = 0, i.e. 

a 

n > , (3.26) 


Following Eq. fl2.10p . in general we dehne 


{0\T^^{q)\0)ren ■= RP (0|T, 

w=0 


fll/ 


(q)|o) 


(3.27) 


For n as in Eq. fl3.26p , consider the expression on the right-hand side of Eq. 03.241) ; 
for any even spatial dimension d this expression is regular in u = 0, so that we 
can apply the zeta regularization in its restricted version to obtain the renormalized 
stress-energy VEV. On the other hand, for odd d we must resort to the extended 
version of the zeta technique, since the function under analysis has a simple pole in 
u = 0. 

Because of the pole singularity, the procedure of evaluating the regular part of Eq. 
03.24P in M = 0 in the case of odd d implies the appearance of a logarithmic term in 
T in the integrand (0). Simple but rather lenghty computations give the following 

^As a matter of fact, for computational simplicity, we will always choose the smaller n G 
{0,1, 2,...} fulfilling Eq. (I3.26p . 

^In fact 

.^u/2 _ gf Inr _ 1 ^ In.,. _(_ (^(y2) for U — >• 0 . 
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results, for d either odd or even and /i, G {0, r,9i}: 


(0|7;,(q)|0),en = + M,,fcT«(r)) , where (3.28) 



+ lnrpW(r;r)J , 

Mi^^k ■= 7sM + 21n(— I 


{'Jem — 0.577 is the Euler-Mascheroni constant). In the above vju){T-,r) and 
Vj^J{T]r) are suitable functions determined by 'H|i“^(r;q); these functions are in 
fact polynomials in of the form 


n+1 

= (ae{0,l}), (3.29) 

i=0 

for some smooth functions '■ [0,+cx)) —R (i G {0 ,...,?t. + 1}). Let us 

stress that 

(^ ; ’^) = 0 , T’Pjir) =0 for d even , (3.30) 

a fact corresponding to the previous comments on the logarithmic terms. 

Next note that, in consequence of the remarks of subsection 13.41 the renormalized 
VEV (0|7)jy(q)|0)ren only depends on the parameter k through the coefficients 
and in the first equation of (13.281) . In particular, the functions T^“^(r) (a G 
{0,1}) introduced therein do not depend on k and we can evaluate them computing 
the integrals in Eq. fl3.28p numerically, for any fixed r G (0, +cxd). 

To conclude, following the considerations of subsection 12.31 (j§), we define the con¬ 
formal and non-conformal parts of the functions r TjiJ{r) respectively 
as 


rn{a,0) 

^ flU 



(aG{0,l}) (3.31) 


The logarithmic term proportional to u is not relevant in the case of even d, where analytic 
continuation exists up to u = 0 and it is simply obtained setting it = 0 in (|3.24p . On the contrary, 
for odd d we must take the regular part at u = 0 of the expression (|3.24l) and the above term 
Inr contributes to it, since it is multiplied by a term proportional to 1/u that comes from the 
M —^ 0 expansion of the factor l/( “~^~^ + l)...( ’^~^~^ +n) in (13.241) . Also recall the statement in 
point i) of subsection 13.41 

^In particular, we use the same convention as in Part I (see Eq. (2.24) therein or Eq. (12.131) in 
the present paper): 

(> = conformal , I = non-conformal . 

Compare Eq. (13.311) with Eq.s (2.25) (2.26) in Part I. 
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is the conformal parameter, see Eq. fl2.14l) ). In snbsections 14.111475] we present 
the fnnctions V^J {a G {0,1}) and the graphs (obtained via nnmerical integration) 
for each one of the functions in Eq. fl3.3ip . for d G {1, 2, 3} respectively. 

3.6 Asymptotics for r = A:|x| —>• 0^". Hereafter we present a method which we 
will use to determine the asymptotic expansion for the renormalized stress-energy 
VEV in the limit r = /c|x| —)■ O'*". 

Let us consider a function F : (0, -|-cx)) —)■ R dehned by 

r+oo 

F(r) := / dr , (3.32) 

Jo 

where h(r) is a positive bounded function, while P{T-,r) is a polynomial of degree 
A in r^; more precisely, we assume that 

N 

Pir^r) = '^Pi{T)r^\ (3.33) 

i=0 

for some integrable functions Pi : (0, -|-cxo) —j- R {i G {0,..., A}). We claim that the 
Taylor expansion of the exponential in Eq. fl3.32p implies an expansion 

N 

F{r) = ^ tti d- PAr+i(r) , PAr+i(r) = for r —)■ 0+ , (3.34) 

1=0 

where the coefficients Oj G R are dehned by 

ai ■= (■_ PI / Pji'r) (^ e {0,..., A}) ; (3.35) 

j =0 d)- Jo 

more precisely, the remainder term is such that 

|PAr+i(r)| ^ for all r G (0,-|-cx)) , (3.36) 

^ -I ^ + 00 

Cn+1 := . (3.37) 

See Appendix for the derivation of the above results. 

Consider the expression fl3.28p for the independent components of the renormalized 
stress-energy VEV; recalling that V^J and vj^J are both polynomials in of degree 
n-|-l (see Eq. fl3.29p ). we note that each component (hxed p, z/) has the form fl3.32p 
fl3.33p with 


/z(r) = tanh(r) , 

P(r;r) = V) + (^K,fc + lnr)Pyj(r;r)j , A = n + 1. 


(3.38) 
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Thus, using the above approach, we can infer the expansion for the renormalized 
stress-energy VEV in the limit r = k\x.\ —)■ O’*"; we defer the explicit results for 
d G {1, 2, 3} to SectionlU where the integrals in Eq. fl3.35p are evaluated numerically 
for the choices of h and P under analysis. 

3.7 Asymptotics for r = fc|x| —>■ -|-oo. Let us discuss some general methods 
allowing to determine the asymptotic behaviour for certain types of integrals; these 
methods will be employed to determine the asymptotic expansion of renormalized 
stress-energy VEV in the limit r = fc|x| —)■ -foo. 

Consider a function F : (0, -|-oo) —)■ R defined by 


F{r) := dv e Q{v,r) (A > —1) 

Jo 

where Q{v,r) is a polynomial in of the form 

Q{v, r) = r) -I- r) Inn , 

N 

(a G {0,1}) , 


(3.39) 


(3.40) 


i=0 


for some smooth integrable functions : [0,1) —)■ R (f G {0, ...,A}). The 

asymptotic expansion of F{r) for large r can be obtained using the standard theory 
of Laplace integrals [71II31EI1I22]; in the present subsection we only report the main 
results of this analysis, making reference to Appendix [A] for the proofs. 

For any K G {1, 2, 3,...} and any vq G (0,1), we have: 


N K+i-l 


A(r) = ^ X] ^A,m{vo,r) + Bi^rn{vo,r)\n.r‘^ 

i=0 m=0 

RK+iivo,^) = Inr'^) 




for r —)• -I-CX3 . 


(3.41) 


In the above expression, the coefficients Ai^rn{vo,r), Bi^rn{'Vo,r) (for i G {0, ...,A}, 
m G {0,1, 2,...}) are defined by 


Ai^rn{vo,r) := ql%'y{m+X + l,Vor^)+ q'i"^Q{m + X + l,Vor^) , (3.42) 

d^ 




Bi^rn{vo,r) := -qlllj{m + X + l,Vor^) , q^l : = 




q^\v) (aG{0,l}), 


D = 0 


where 7( , ) denotes the lower incomplete gamma function and we put 

g(s,^):= [ dw e~'^w^~^\nw for all 2 ;G (0, -|-cxd), sGC with 3fJs > 0 . (3.43) 
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Concerning the remainder term RK+i{vQ,r)^ we can give qnantitative estimates of 
the form 

\RK+i{yQ,r)\ ^ {Fk +1 + Gk+i Inr^) for all r G (0, +cx)) ; (3.44) 


these involve some constants Fk+i-,Gk+i > 0, that can be determined explicitly 
following the indications of Appendix lAl 

Moreover, from the asymptotic behaviour of the incomplete gamma 7 ( , ) and of 
the function g( , ), we infer 


N K+i-1 




2 = 0 771=0 

SK+i{vo,r) = for r^+oo , 


(3.45) 


where we have defined the real coefficients (for iG{0, mG{0,1,2,...}) 


i^rn ■= r(m+A + l) + 'ijj{m+X + l) g, 

A,m := -r(m+A + l) ; 


R) 


(3.46) 


in the above definitions, 'ip{s) := 9slnr(s) denotes the digamma function and the 
qi^miQilm S'l’e as in Eq. fl3.42p . In principle, the remainder term SK+i{vo.,r) in Eq. 
fl3.45p could be evaluated quantitatively combining Eq. fl3.44p with some known 
estimates about the incomplete gamma functions; the related computations would 
be very tedious, and will be avoided here. 

Let us now consider expression (I3.28p for the renormalized VEV of the stress-energy 
tensor and perform inside the integrals appearing therein the change of variable 


T = arcth (n) , n G (0,1) 


(3.47) 


(“arcth” denotes the hyperbolic arctangent); in this way we obtain 



dv 


(arcth n) 


d+3 


(arcth v ; r) -|-ln(arcth v) (arcth v ; r) 




^ ^ Pil^Harcthn ;r) . (3.48) 

Since vjn) and V^J are polynomials in of degree n-l-l (see Eq. fl3.29p ). we can 
re-express each component of the renormalized stress-energy VEV (fixed /x, z/) in the 
form (13.3911 ■ where Q{v,r) is as in Eq. fl3.40p with 


A := n — 


d + 3 
2 


N = n + 1 ] 


(3.49) 
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(3.50) 


1 

l — v‘^ 


g(0)(^.^) _ 


/ arcthf 

V V 


d+3 


Vj^J (arcth v ; r) 


ln|^ 


/ arcth V 


+ M^^k]vl^J{aTcthv;r) 


Q^^\v;r) := arcth (arcthv ; r) . (3.51) 


In this way we can derive asymptotic expressions analogous to fl3.41l) (and fl3.45l) i 
for the renormalized stress-energy VEV in the limit r = fc|x| —>■ -l-cxo; we report in 
Section m the explicit results of this computation for the cases with d G {1,2,3}, 
respectively. 


3.8 The total energy. Let us recall that the total energy can always be expressed 
as the sum of a bulk and a boundary contribution as in Eq. fl2.36p : in the following 
we are going to discuss these two contributions separately. 

Let us hrst consider the regularized bulk energy; according to Eq. fl2.37p . this is 


E^ = - Tr ^ V . 

2 

This is connected through Eq. (I2.23p to the heat trace := Tre“^ that, ac¬ 
cording to Eq. (I3.9p . has the form 

m = with (3.52) 


it is patent that the map 1 1 —)■ H{i) is smooth on [0, -|-cxd) and exponentially vanishing 
for I —)■ -l-cxD. Then, following the general considerations of subsection 12.81 (see, in 
particular, Eq. fl2.28p for the trace Tr^“*), we obtain for the regularized bulk 
energy 


= 


-ir 


^+oo 




di i 




(3.63) 


The above relation holds for any n G (1, 2, 3,...} and the integral appearing therein 
converges for any complex u with 


IRu > 2{d — n) -|- 1 ; (3.54) 

thus, for any integer n > d + 1/2, Eq. fl3.53p gives the analytic continuation of E'^ 
in a neighborhood of m = 0. Since here no singularity appears, we can obtain the 
renormalized bulk energy simply by setting m = 0 in Eq. fl3.53p : making again the 
change of integration variable r := k'^t (see Eq. fl3.15p ). we infer 
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E' 


^ren 



(3.55) 


The above expression is fully explicit and holds for any spatial dimention d; in the 
subsequent Section 0] we are going to evaluate numerically the integral appearing in 
Eq. fl3.55p for d G {1, 2, 3}, making the minimal choice n = d + 1. 

Let us point out that in the case under analysis, where the potential is assumed to 
be isotropic, we could also derive an alternative representation for the bulk energy 
Eren terms of the (analytically continued) Riemann zeta function; we discuss this 
topic in detail in Appendix |Bl Here we prefer to focus the attention on the approach 
fl3.55p for the computation of since it is more general; in fact, it could be 
employed straightforwardly also when the background potential is not isotropic (j§). 
As pointed out in the Introduction, the computation of the renormalzied bulk energy 
was also performed by Actor and Bender in [2]; these authors employed a global zeta 
approach based on ad hoc results on the analytic continuation of some particular 
kind of “generalized” zeta functions (discussed in detail in PQ), befitting the con¬ 
figuration under analysis. Let us stress once more that here, instead, we derive the 
required analytic continuation applying a general procedure in a mechanical way. 
Our method and the one of |2] could be proved to be equivalent in any dimension, 
yet the general discussion is too involved to be reported here; in subsection 14.31 we 
check by direct comparison that, in the case d = 3 (Q), the numerical result obtained 
using prescription fl3.55p agrees with the one derived in [2]. 

Now, let us move on to discuss the boundary energy. In all the cases treated in 
Parts I and II this boundary contribution was easily found to vanish identically in 
consequence of the boundary conditions. Here the situation is not so simple: since 
the spatial domain is unbounded (indeed, H = R'^), we must resort to the procedure 
pointed out in subsection 12.101 and intend the integral over dfl in dehnition fl2.38p as 
the limit of integrals over the boundary of suitable, bounded subdomains {^e)e=o,i,2,... 
(see the comments below Eq. fl2.38p L Due to the spherical symmetry and to the 
characteristic scale of the problem under analysis, it is natural to choose for VLi the 
balls with center in the origin and radius ^/k^. 


:= R^/]^(0) = {x G I |x| ^ for £ G {0,1, 2,...} . (3.56) 


Then, Eq. fl2.38p for the boundary energy reads 



(3.57) 


®Let us recall that in the present paper we have chosen to omit the analysis of cases with 


unisotropic background potentials only for simplicity. 

^As a matter of fact, this is the only model discussed both here and in [5]. 
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where indicates the {d — l)-dimensional spherical hypersurface, 

while c?ry denotes the derivative in the radial direction with respect to the variable 
y. In terms of the set of rescaled spherical coordinates fid.lip . we have 


a i\:(t;x,y) 


y=x 


A/27rsinh(2fc2t) 


{k r) tanh(fc2t) e”"'; (3.58) 


then, using the integral representation fl2.22p for the Dirichlet kernel in terms of the 
heat kernel (again, with the change of variable r := and noting that the above 
expression for the latter does not depend on the angular variables (^i, ...,6*d_i) (j^), 
we infer 


Bu 


k 

2^-'r(f)r( 


u+l \ 
2 > 



(3.59) 


where, for brevity, we have put 


r+oo 

03“:= lim / drr"^ 

£^ + 00 Jq 


tanh r 


^d/2 itanhr\ 
(sinh(2r))'^/^ / 


(3.60) 


For any d G {1,2,3,...}, £ ^ 0, the expression within the round brackets in the 
above equation is an analytic function of r on [0, +cxd) and vanishes exponentially 
for T —)■ +CXD. In consequence of this, we easily infer that the integral in Eq. fl3.60p 
converges for any complex u with 


^u>d-3 . (3.61) 

Within this region, by Lebesgue’s dominated convergence theorem (Q) we can take 
the limit under the integral sign in Eq. fl3.60p . and the conclusion is 

03“ = 0 for sRm > d - 3 . (3.62) 

®We also use the fact that the area of the (d — l)-dimensional spherical hypersurface of radius 


r is 


Area(^si^ 


fidJY) 


r 


d-l 


^Indeed, notice that, for any f > 0, there holds 


f tanh T 


£d^d/2 i|;a,nhT 


u-d+i ( 1 tanh r 

< r 2 f 


(sinh(2r))^/2 

where the right-hand side is Lebesgue-integrable over (0, -l-oo); besides, 


(sinh(2r))‘^/^ / \2etanh(T) 


hm g-^=tanhr 


^^+oo 


(sinh(2r))^/2 


= 0 for all r G (0, -|-oo) . 
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Summing up, by analytic continuation we conclude that the renormalized boundary 
energy vanishes: 


^ren _ 


= 0. 


ti=0 


(3.63) 


4 The previous results in spatial dimension d G {1, 2, 3} 


4.1 Case d = 1. Let us hrst recall that in this case the conhguration is symmetri¬ 
cal under the reflection —x^; on any one of the intervals (0, -|-cx)) and (—cxo, 0) 

we can use the coordinate q = r := k\x^\ G (0, -|-cx)) (see Eq. fl3.12p L 
After lenghty computations, based on Eq.s fl2.161l2.18|) and on the general consid¬ 
erations of subsection 13.51 we obtain for the regularized VEV of the stress-energy 
tensor an expression of the form fl3.18p . More precisely, we have 


^ n2 / ^+00 

(0|^M^.(q)|0) = (^) / dr (/r, z/G {0, r}), (4.1) 


r(^) _ 

where (recall that dependence on ^ is understood) 

^TO(^;q) := 


Ai(r,r) 


(l-u)(l+40 + (l-40 
2r (1-1-4,^ cosh 2r) 


2r 


sinh 2r 


1 + 


sinh 4r 
2 cosher 


(1-40 1-W + 


2r 


H-cosli r 


(4.2) 


(4.3) 


H'“>(T;q):=yl,(r,r) 

and the off-diagonal components hIm) (r ; q) (/i 7^ i^) vanish identically; in the above, 
for simplicity of notation we have set 


Ai(r,r) := 


tanh r 


2t 


16-0^ 


sinh(2r) 


(4.4) 


The above expressions for the components of the tensor (r ; q) are easily seen to 
possess the features anticipated in Eqs. fl3.22p fl3.23p and in the related comments. 
Thus, according to the general framework developed in subsection 13.51 we can obtain 
the analytic continuation of (0|T“^,(q)|0) given in Eq. fl4.ip integrating by parts n 
times the integral therein, with n > 1 (see Eq. fl3.261) 1. Let us fix 


n = 2 , 


so that Eq. fl3.24p gives 
(0|T;,(q)|0) 


e 1 

r(!!±i) (f-i)f 



+00 


dTri5;wW(T;q). 


(4.5) 
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(4.6) 
































It is evident that (0|T“j^(q)|0) has a simple pole at m = 0 (this reflects a general 
feature of the cases with odd spatial dimension, already indicated in subsection 13.51) . 
Thus, the renormalized stress-energy VEV is given by the regular part in u = 0 of 
the function in Eq. fl4.6p (see Eq. fl3.27l) h with some effort, we can re-express this 
quantity as in Eq. 03.281) . i.e. 

(0|f,„(q)|0)„„ = ^(rfyr) + AUkTli}(r)) , (4.7) 




{r) := dT e-’-’ *“>'" P<“„> (t ; r) + In T (r ; r) 


^+00 


T^l\r):= I dre 


—r^ tanh r 


where 




^ tanh i 


TT 


'PlLHT\r) , := 7B„+21n('y V 


3?Wi“>(T;r) + 2a„ S|Wi:>(r;r) 


u=0 


fii/ 


P<pT;r) ^ e' “>■" a?WS(r;r) 


(4.8) 


TT 


It is readily found that Vj^J, are polynomials of degree iV = 3 in r^. We can 
evaluate numerically the integrals in Eq. (14.7p . distinguishing between the conformal 
and non-conformal parts (>, B of each component: see Eq. 03.3ip . recalling that for 
d = 1 we have (see Eq. 02.14p l 

= 0 . (4.9) 

The forthcoming Fig.s [IHshow the graphs of the functions 


r I-)- T^“’^^(r) for/i = z/G {0,r}, a e {0,1}, X e {0, • (4-10) 



Figure 1: d = 1: graphs of Tqq’*^^ and Tqq’^^ . 
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Figure 2: d = 1: graphs of Tqo’"^ ^^^1 Tqq’"^ • 



Figure 3: d = 1: graphs of Trr’^'^ and Trr’^'^ 



Figure 4: d = 1: graphs of Trr’"^ and Trr’"^ • 
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Let us pass to evaluate the small and large r asymptotics of the functions in 
Eq. fl4.10p . On the one hand, using Eq.s (13.32113.35]) (with iV = 3) we obtain, 
for r = k\x^\ —)■ O’*", 


= -0.0153 + 0.0164- 0.0796 + 0.0262 r® + 0(r®) , 

= 0.0398 r2 + 0(r«), 

T(S°’")(r) = 0.2121 - 0.3766 + 0.2356 - 0.0903 r® + 0(r®) , 

hl'V) = 0(r*) i 

= -0.0153 - 0.0164+ 0.0265 - 0.0052 r® + 0(r®) , 

rji’0)(r) = -0.0398 + 0(r®) , 

(4 12) 

T^(0-")(r) = -0.0002 - 0.0001 + 0(r®) , 

r»’')(r) = 0(r*) , 

The numerical coefficients appearing here and in other small r expansions are ob¬ 
tained from numerical computation of the integrals in fl3.35p . 

On the other hand, using Eq.s fl3.39p (13.451) (13.461) (with 77 = 3, A = 0) we obtain 
the following asymptotic expansions, for r = k\x^\ -|-oo: 


p(o,0) 

-00 

p(iiO) I 
-00 ' 

40, 


(r) = -^(1„+ 7EM -H) ^ Inr^) . 


r) = -h 0{r ® In r^) , 

Stt 


(4.13) 


Tr\r) = 

41,1 


27rr^ 37rr® 

2 


-|- 0{r ® In r^) , 


T^^'~\r)=0{r-^\nr^) ; 


TS-0\r) = ^ (In + 7 EM - l) + ^ + + 0(r-« In r^) . 


7 


= - -hO(r ®lnr^) , 

87r 

T^^°’"^(r) = 0(r“®Inr^) , 
Trl'*\r) = 0(r“®lnr^) . 


(4.14) 


Let us now discuss the bulk energy; using the expression (I3.55P with n = 2 and 
evaluating numerically the integral appearing therein, we infer 


= (0.0430546469 ± 10"^°) k . 


(4.15) 
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4.2 Case d = 2. Consider the general framework of snbsection 13.51 in this case 
we nse the (rescaled) polar coordinates q = (r, 6*i) G (0, +cx)) x [0, 27r), fnlhlling (see 

Eq. dMlD) ( 0 ) 

kx^ = rcos9i, kx‘^=rsm9i. (4.16) 

Proceeding as in the one dimensional case, with some effort we obtain the following 
integral representation for the zeta-regnlarized stress-energy VEV (compare with 
Eq. fl3.18p and recall that dependence on ^ is nnderstood): 


(0|E;).(q)|0) = (k) dr r“ 2 +^?fj,“j(r;q) (/i, i/G {0, r, di}) , (4.17) 

where the tensor is diagonal and, concerning the diagonal components, we have 

: = 


dl2(r,r) 


-(l-w)(l+40 + (l-40 


2r 


sinh 2r 


2 + C 


:= A2{T,r) 


se 


T 




tanhr 

T 


(1 — 4^) ( 1—u + r 


sinh 4r 
2 cosher 

2 2r 
cosher 


8^' 


tanhr 

In the above, for simplicity of notation we have pnt 


(1 — 4^) il — u -f r' 


2r cosh2r 
cosher 


A 2 {r,r) := e 
64 TT 


—r^ tanh r 


2r 


sinh 2r 


(4.18) 

(4.19) 

(4.20) 

(4.21) 


Again, the featnres indicated by Eq.s fl3.22p fl3.23p (and related comments) are all 
possessed by the expressions (14.18114. 2 II) for So, we can analytically continne 

in u the expression in Eq. (I4.17p integrating by parts n times, with n > 3/2 (see 
Eq.s fl3.24p fl3.26p h we choose 

n = 2 , (4.22) 


giving 


(0|T"(q)|0) = 


k^ 






+ OD 

drrt-ia2-HW(r;q) . (4.23) 


Patently, each component of (0|T/I|0) is regnlar at m = 0; thns, we can obtain 
the renormalized version of the stress-energy VEV by simply pntting m = 0 in Eq. 

course, the spatial line element in this coordinate system reads = k~'^{dr"^ +r^dO^); this 

determines the Christoffel symbols in Eq. (12.421) for the derivatives Dij. 
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fl4.23p (we already noticed this property in general for the cases with even spatial 
dimension d; see below Eq. fl3.26p h In conclusion, 


(0|T;„(q)|0)„„ = fc‘’r<“>(r) , 


(4.24) 



—tanh r 



r; r 


), 




tanh r 



(compare with Eq.s fl3.281l3.30p h also in this case, we easily check that V^J is a 
polynomial of degree = 3 in r^. Next, we proceed as we did in the previous 
subsection for the case of spatial dimension d = 1: we evaluate numerically the 
integrals in Eq. fl4.24p and separate the conformal and non-conformal parts 0, B of 
each component (see Eq. fl3.3ip h noting that Eq. fl2.14p gives 


6 = 


1 

8 ■ 


(4.25) 


The forthcoming Fig.s[MZlshow the graphs of the functions 

r ^ T^°’X)(r) (/i = z/G{0,r}) , {k/rfT^X^{r) for x G {0,"} • (4.26) 
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7'(0.O) ^ 00 

^00 



Figure 5: d = 2\ graphs of Tgo’*^^ and Tqo’"^ • 



Figure 6: d = 2\ graphs of Trr’^'^ and T^r’"^ • 



Figure 7: d = 2: graphs of {k/r^Tg^f^'^ and {k/r)‘^Tg^^^. 


28 




















Now, let us consider the small and large r asymptotics of the functions in fl4.26p . 
On the one hand, Eq.s (13.32113.35|) (with N = 3) give, for r = k\x.\ —)■ 0+, 

= -0.0017 - 0.0134 - 0.0154 + 0.0027r® + 0(r®) , 

T(SJ’")(r) = 0.1649 - 0.1069 r2 +0.0516 r^-0.0141 r® + 0(r®) ; 


= -0.0010 + 0.0207+ 0.0114- 0.0013 r® + 0(r®) , 
TjO’")(r) = -0.0806 + 0.0267- 0.0133 + 0.0018 r® + 0(r®) ; 

{k/r^T^^^gfir) = -0.0010 + 0.0140 + 0.0153 - 0.0027r® + 0(r®) , 
{k/rYT^^^g"\r) = -0.0806 + 0.0802 - 0.0440 + 0.0124r® + 0(r®) . 


(4.28) 

(4.29) 


Again, the above coefficients were obtained calculating numerically the integrals in 
Eq. dMSD. 

On the other hand, using Eq.s (13.391) (I3.45P (I3.46P (with K = 5, X = —1/2; note that 
no logarithmic term Inr^ appears in this case) we obtain the following asymptotic 
expansions for r = /c|x| +oo: 




19 


+ 0 (r 


-9^ 


127r 25607rr® 

^ ^ +0(r-9); 


^orW = 7r:: + 


47rr 327rr® 


T^r'^\r) = ^ + 


17 


127r 487 rr 25607rr® 

TS'^ir) = ~ + -At + O(r-0) ; 


+ 0(r 


-9^ 


Wr)^T^°£\r) = 


Anr 3271 
1 


+ 


33 


127r 967rr ' ' 25607rr® 


+ 0(r 


- 9 \ 


(k/rynTi^) = 


377 


+ o{r-y . 


(4.30) 


(4.31) 


(4.32) 


As for the bulk energy, using again the expression (I3.55P with n = 3 and evaluating 
numerically the integral appearing therein, we obtain 


Even ^ _(o.oi80207591 ± 10“^®) k . (4.33) 


4.3 Case d = 3. In this case we use the coordinates q = (r, ^i,^) G (0,+cxo) x 
( 0 , 7 r) X [0,27r), which are related to the Cartesian coordinates x = (xi, 2 : 2 ,Xs) via 
(see Eq. (13. lip ) (0) 

kx^ = rcos9i, fc = r sin 6*1 cos 6*2 , fc = r sin 6*1 sin 6*2 . (4.34) 

this case the framework of subsection 12.111 must be employed using the spatial line element 
di^ = k~^{dr^+r^{d6l+sm^6id92)) and the corresponding Christoffel symbols. 
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Similarly to what we did in the previous two subsections, after lenghty computations, 
we can express the zeta-regularized stress-energy VEV as (compare with Eq. flS.lSp 
and recall that dependence on ^ is understood) 


^ ^4 / \„ r+oo 

(0|^M:.(q)|0) = (fc) dr (/r, i/G {0, r, 0i, ^a}) , (4.35) 


where 'H'ftv is diagonal and we only have to consider the independent components 




-(l-^i)(l+40 + (l-40 


2 r 


sinh 2r 


3 + r" 


sinh 3r — sinh r 
cosh T 


'H^)!^(r;q) := 


413(7-, r) 




r 


tanhr 


(1-40 1-W + 


2 r 


sinh 2r 


l-|-2r^ tanhr 


r V 


k) 


8 ? 


r 


tanhr 


(l-40fl-« + (- 


2 r 


cosh 2r 


(4.36) 


(4.37) 


(4.38) 


^sinh2ry 

(for the remaining diagonal component, i.e. , see Eq. fl3.22p h In the 

above, for simplicity of notation we have put 




2 r 


sinh 2r 


3/2 


(4.39) 


Also this time, the expressions fl4.36114.39P for possess the properties indicated 
in Eq.s fl3.22p fl3.23p (and in the related comments); thus, we can obtain the analytic 
continuation in u of the expression in Eq. fl4.35p integrating by parts n times, for 
any n > 2 (see Eq. 03.261) h For dehniteness, we £x 


n = 3 , 


(4.40) 


so that Eq. 03.24p reads 


( 0 |f;„(q)| 0 ) 


k* 1 

r(!!±i) (f- 2 )(| 




+CXD 

<;TTia?WM(r;q). 


(4.41) 


As in all cases with odd spatial dimension, the analytic continuation of the reg¬ 
ularized stress-energy VEV given in Eq. 04.4ip has a simple pole in m = 0 (recall 
subsection 13.5p . In consequence of this, we have to adopt the extended version of the 
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zeta approach to define the renormalized VEV (0|T^i,(q)|0)ren, taking the regular 
part in M = 0 of Eq. fl4.4ip (see Eq. fl3.27p h with some effort, we obtain 


(0|f„„(q)|0)„„ = *:*(r<?(r) + A4,»rW(r)) , (4.42) 


r'?(r):=/ T’(»)(T;r) + lnTpW(r;r) , 

r™(r) := r°°dT P'pT; r) , A4., := 7 em +2 \n(^'] 


where 


Pi°’(r;r) := 


tanh r 


40r 


3 q) + 4 (t; Q) 


4Z=0 


e’'’*""' a?W<2(T;q) ; 


(4.43) 


this time, we readily infer that Vj^), are polynomials of degree iV = 4 in r^. 
Now, we evaluate numerically the integrals in Eq. (I4.42p and distinguish between 
the conformal and non-conformal parts 0, B of each component; once more we refer 
to Eq. (13.3ip . recalling that for d = 3 we have (see Eq. (I2.14p ) 


^3 — c • 

D 


(4.44) 


The forthcoming Fig.s ISlfT^ show the graphs of the functions 


r^T(“’")(r), T(“’X)(r), forae{0,l}, xe{0,B} • (4.45) 


-,(a,x) 
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Figure 8: d = 3: graphs of Tqo”^^ and Tqq’^^ 


0.08 



r 


0.06 

0.04 

0.02 


Figure 9: d = 3: graphs of Tqo’"^ ^oo’"^ • 


10 


r 


J’(hO) 



Figure 10: d = 3: graphs of Trr’^'’ and Trr’^^ 
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Figure 11: d = 3: graphs of Trr’*^ and • 



Figure 12: d = 3: graphs of and . 



Figure 13: d = 3: graphs of {k/r)‘^TQ^^^ and {k/r)‘^TQ^^^ . 
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In conclusion, let us discuss the small and large r asymptotics of the functions in 
fl4.45p . On the one hand, Eq.s fl3.321l3.35|) (with = 4) yield, for r = fc|x| —)■ 0+, 


To^S’^)(r) = -0.0047 - 0.00240.0028 r^+ 0.0006 r®- 0.0001 r®+0(r^°), 
T^^’^\r)= -0.0016r^+0(r^°) , 

T(S°’")(r)= -0.0143 - 0.0468r2+ 0.0134r^- 0.0033r®+ 0.0007r® + 0(r^°), 
^oo’"^W= 0.0380+ 0(ri°) ; 


= -0.0016 + 0.0039r2- 0.0003/- 0.0005r®+0(r^°) , 

Tji’0)(r) = 0.0016/ + 0(r^°) , 

= 0.0095 +0.0188/-0.0038/+0.0007/-0.000l/+0(r^°), 
Tji«(r) = -0.0253+ 0(ri°) ; 

(k/rf (r) = -0.0016 + 0.0023 + 0.0004 / - 0.0006 / + 0.0001 / + 0(r^°) , 

(k/rfT^^^;)^\r) = 0.0016/ +0(r^°) , 

(fc/r)2Tj°*(r) = 0.0095 + 0.0375/ - 0.0115/ + 0.0030/ - 0.0006/ + 0(r^°) , 
(k/rfTg^^0"\r) = -0.0253 + 0(r^°) . 

(4.48) 

On the other hand, Eq.s fl3.39p fl3.45p fl3.46p (with 77 = 4, A = 0) allow us to infer 
the following asymptotic expansion, for r = fc|x| —)■ +cx3: 


rj-i( 0 , 0 ) 

-^00 


(r) 


rp(l, 0 ) 

-‘■00 


(r) 


rp(o,m) 

-‘oo 


(r) 


rpilM) 

-‘oo 


(r) 


-/ (lnr^-^7EM+ 1) + /s;! + 0(r-»lur=) , 
+ + 0(r~^ In rO ; 


(4.49) 


ry»>(r) 

Ti'/\r) 


^ + 0(r-«lnr2) , 

+ (In r'n7EM) + + 0(r-« In r^) , 

-/ + 0(r-''lnr") ; 


(4.50) 
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+ 0(r ®lnr^) , 


(4.51) 


Finally, Eq. fl3.55p with n = 4 and numerical evaluation of the corresponding integral 
allow us to derive the bulk energy 


Even ^ -(0.0078607119 ± 10"^°) k . 


(4.52) 


This result agrees with the one obtained by Actor and Bender [2] using a different 
method ff^b 

Acknowledgments. This work was partly supported by INdAM, INFN and by 
MIUR, PRIN 2010 Research Project “Geometric and analytic theory of Hamiltonian 
systems in hnite and inhnite dimensions”. 


check this, one must compare the numerical value reported in the above Eq. (14.521) with 
the one reported in Eq. (4.4) of [2]. Let us stress that conventions different from ours are used 
therein. In fact, using our language, the bulk energy is formally defined in [2] as E := 
while our general prescription (12.3711 gives E = ^ moreover the parameter a of [2] and our 

parameter k are related by a = '/2 k. Summing up, the “total energy” derived in [5] has to be 
multiplied by 1/2 in order to obtain our 
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A Appendix. Asymptotic expansions for certain 
integrals 

Hereafter we show how to obtain the expansions (13.341) (13.411) . holding respectively 
for the functions dehned via the integral representations fl3.32p fl3.39p . The proofs 
are given in subsections lA.ll and IA.31 subsection IA.2l is an interlude on gamma-type 
functions, useful in view of subsection IA.3I 

A.l Derivation of Eq.s (I3.34H3.5711 . Let us consider the framework of subsec¬ 
tion [3Tl where 

^+00 ^ 

F{r) := / dr P{t ;r) , P{t ;r) ='S^Pii'^) ^ (0,+ 00 )) , (A.l) 

>^0 *=o 

{h a positive bounded function, pi {i G {0, some integrable functions). 

To evaluate P{r) we start from Taylor’s formula with Lagrange remainder for the 
exponential; this ensures that, for any z G (0, -fcx)) and any M G {0,1, 2,...}, 


M 


e-" = 




Pm+i{z*) 


m\ 

m=0 

1 


(M+l)! 


- 1 ) 


M+l 


(A.2) 


2 = 2* 


(M+l)! 




for some + = z^{z) G (0,2;); in particular, we have 

Ipm+i(+)| ^ ^ 


(M+l)! ■ 


(A.3) 


Substituting expansion (1A.211 with z = r^h(r) and M = N — i into Eq. (lA.ip . we 
readily infer 


N 


F{r) = 


i=0 


'N-i 

E 

m=0 


-1) 


m r+co 


m\ 


dTpi{T) (h(r))"* + 


b+oo 


dr pi(r) (h(r))^ pN-i+i{z^{r,r)) 


(A.4) 


with +(r, r) G (0, r^h{T)). Introducing a new summation index i such that m = j—i, 
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and performing the exchange i j (jlj), the above relation yields 


N 


F{r) = aj + Rn+i jr) for r e ( 0 ,+cx)) , 


(A.5) 


i=0 


where the coefficients Oj G R are dehned as in Eq. fl3.35p . while the remainder term 
RN+i{r) is 


RN+i(r) := 


/ ^ /• + 00 

Vi=o -^0 


drpiir) (h(r))^ *+ViV-i+i( 2 :*('r, r)) . (A. 6 ) 


Using the estimate flA.3p . we readily infer the nniform bonnd, holding for all r G 
( 0 ,+oo), 


RAr+i(r)| ^ 


N 

E 


^+00 


^ (iV-z + 1)! 7o 


dr \pi{T)\ (h(r)) 


^^-*+1 i ^ 2 ( Af + l ) 


(A.7) 


. i =0 


this proves Eq. (I3.36p . with the expression (13.371) for the constant Cn+i therein. 

A.2 The incomplete gamma functions and the integral g(s,x); asymp¬ 
totics and bounds. Let us hrst consider the “lower” and “upper” incomplete 
gamma functions, respectively dehned as 

7 ( 5 , z) := f dw e~'^ for all s G C with 3fJs > 0, z G (0, + 00 ) , (A. 8 ) 


^ + 00 


r(s, z) := J dw e for all s G C, 2 ; G (0, +cxd) . (A. 9) 

Hereafter we report some well-known properties of these functions (see [20] , Chapter 
8 ), to be used in the following subsection. First of all, there holds 


7 (s, z) + r(s, z) = r(s) (s i { 0 , - 1 , - 2 ,...}), 


(A. 10) 


where r( ) is the Euler gamma function. Concerning the lower incomplete gamma, 
there hold the relations 


7 (s-|-l, z) = 57 ( 5 , z) — e , 7 ( 5 , z) = r(s) + 0{z °°) for 2 : +cx) ; 

0 ^ 7 ( 5 , z) ^ r(s) for s G (0, + 00 ) . 


(A.ll) 


^^Notice as well that, for any family it is 

N N N i 

j—O i—j j—Q 
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On the other hand, the upper incomplete gamma fulhlls 

r(s + l, z) = s r(s, z) + , r(s, z) = 0{z~°°) for z —)■ +oo ; 

0 ^ r(s, ^ r( 2 ;) for s G (0, +oo) . 


(A. 12) 


(both here and in Eq. flA.lip . the remainder 0{z °°) indicates a quantity which is 
0{z-^) for all N G {0,1,2,...}). 

Let us move on to discuss the properties of the function 0 (s,a;) dehned via the 
integral representation (see Eq. fl3.43|) ) 

0 (s, z) := / dw Intc for all s G C with > 0, 2 ; G (0, + 00 ). 

Jo 

First of all, since ds{w^~^) = Into, from the above definition and Eq. (lA.Sp it 
trivially follows that 

0 ( 5 , 2 ;) = 7 ( 5 , 2 :) ; (A.13) 

thus, using the recursive relation in Eq. flA.llj) . we can easily infer 

0(5 + 1 , ;2) = 50(5, 2 ;) + 7(5, ;2) — 2 ;® ln ;2 ( 3 fJ 5 > 0) . (A.14) 

Now, let us show that 

0(1,2;) =-e“Mn2; - 7 £;m - r(0,2;) (A.15) 

ilEM denotes the Euler-Mascheroni constant). To this purpose, let us write 

0 ( 1 , 2 ;)= lim 05 ( 1 , 2 ;), 05 ( 1 , 12 ) := [ dwe~'^\nw (5 > 0 ) . (A.16) 

< 5 ^ 0 + Jg 

Consider the regularized function 05 ( 1 , z) and integrate by parts to obtain 


05 ( 1 , 2 ;) = —e * Inx + e In 5 + / dw e '^w ^ ; 


(A.17) 


now, note that (for 0 < 5 < 2 :) 

/*z /*+oo 


^+00 


dwe '^w dw e '^w dw e '^w ^ = r(0, 5) — r(0, 2 :) . (A.18) 

Js Js Jz 

Since r(0, 5) = — In 5 — 'yEM + 0(<5) for <5 —)■ 0+ (see [20], p.l77, Eq. 8.4.15), we have 
05 ( 1 , 2 :) = —Inx — 7 £:m — r(0, z) — (1 — e”*^) ln5 + 0(h) for h —)■ O’*", (A.19) 
which, along with Eq. flA.16p . yields Eq. flA.15D . 
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Let us proceed to prove that there holds the bound 


fl(s,z) ^ -s 


-2 


for all s, z G (0, +oo) ; 


(A. 20 ) 


we are going to give the proof for z < 1 and z ^ 1 separately. On the one hand, for 
z < 1 (and s > 0 ), there holds the following chain of inequalities: 


g{s,z)^ dwe'^w^lnw^ dw \nw = —s , 


(A.21) 


where the first passage follows from the negativity of the integrand (Into < 0 , for 
w G (0,1)), while for the second we used e“'^ ^ 1 for w G (0,1). On the other hand, 
for z > 1 (and s > 0 ), we have 

rl rz rl 

g{s,z)= / dw e~'^w'^~^\nw + / dw Inw ^ / dw \nw (A. 22 ) 

Jo Ji Jo 

since the integral over (1,^;) is positive; then, the thesis flA.20p for 2 ; > 1 follows 
using the same arguments as in Eq. flA.2ip . 

Finally, we prove the asymptotic behaviour 

g(s, z) = r(s) ^jJ{s) + 0{z~°°) for all s G C with > 0, 2 ; —)■ +00 , (A.23) 

where s 1 —)■ 'ip{s) := c?slnr(s) is the digamma function. To this purpose, write 


^+00 


^+00 


0 (s,; 2 )= / dw e ^\nw— dw e ^Intc. 

Jo Jz 

Concerning the hrst integral, we have (see |ll], Eq. 4.352.1) 


r"+oo 


dw e "'w® ^ Intc = r(s) 'ijj{s) {s G C, 3fJs > 0) , 


(A.24) 


(A.25) 


so that Eq. flA.23p follows if we can prove that 


™+oo 


dw e '^w'^ Mntc = 0{z °°) for s G C with 3fJs > 0, z ^ +00 . (A.26) 


Indeed, for any 2 : G (1, +C)o), we have 


0 ^ 


^+00 


dw e ^w^ ^ In w 


^+00 




dw e-^w^^ = T{^s + l,z) , 


(A.27) 


where the second inequality follows from the fact that | Intel ^ w (for w G {z, +cxd) C 
( 1 , +cxd)), while in the last equality we used the dehnition flA.Qp of the upper incom¬ 
plete gamma function (with 3fJs > 0). Now, Eq. (1A.26P follows from Eq. (]A.12p . 
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A.3 Derivation of Eq.s (I3.45|) (I3.46|) . Consider the framework of subsection 
[321 where 


F{r)= I dv e-'^"^v^Q{v;r) (A >-1) , 

Q{v] r) = (n; r) + (n; r) In n , (n; r) = ^ {v) 


N (A.28) 


i=0 


for some smooth integrable functions : [0,1) —?■ R (i G {0,...,A}). Let us 

begin hxing vq G (0,1) and re-expressing F{r) as 


N 


= Y1 r) + 4^^(no, r) + (^^o, r) + (^o, r)^ , (A.29) 

i=0 ^ 

where we set, for i G {0,..., N}, 

pVQ pi 

I-^\vo,r) := / dv v^qf’\v) , I-^\vo,r) := / dv v^qf\v) ] (A.30) 


'VO 


rvQ pL 

J-^\vo,r) := / dn n^g-^^(n) Inn , J-^\vo,r) := / dv v^ql^^\v)\nv . 


'VQ 


Concerning the integrals on (ng, 1), we readily infer the bounds 


4^H^'o,r) 


^ e 




'VQ 


j4^(no,r) 


^ e 


/ dv lii(,j) (u) I 


(A.31) 


' 1^0 


Let us move on to discuss the integrals 4^\ 4^^- For any hxed i G {0,..., N} and 
any M G {0,1,2,...}, consider the Taylor expansions near n = 0 with Lagrange 
remainder 


p.'rLiW - dPSTfi 

Of course. 


M 

m=0 


Aa), 


AO W ^ 0,M+l('^o) (a) 


(A.32) 

for some n* = n*(n) G (0,n) . 


sup |p*,M+i(^*)l ^ TyFTTAT’ 4 m+i(uo) := sup 

ve{o,vo) [M+L)\ ^g(0,„, 


(0,i^o) 


dM+l 




(“) / ^ 

Qi (u) 


(A.33) 
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Now, substitute the expansion flA.321) into Eq. flA.30p and make the change of vari¬ 
able V = wjr^ (with ta G (0, Uo in the (<) integrals appearing therein. Recalling 
the definitions flA.SD and fl3.43p respectively of the lower incomplete gamma 7( , ) 
and of the function g( , ), we obtain (Q) 


M 


Q^l^im+X + l^vor"^) r 2(»n+^+i) + R^Yii^o^r) , (A.34) 


m=0 


jY{vo,r)= (A.35) 

M 

Qim + X + l.VQr"^)--f{m+\ + l,Vor'^)\nr‘^ r) , 


m=0 

where the remainder functions are dehned as 
?{o) 


^U/+i(^o,r) := / dve P*,M+i(n*) 


^lM+i(^'o,r) := / dve ^ 


>0 

"Vo 


( 1 ) 


M+A+1 


(A.36) 


Inn 


Then, using the bound flA.33D and the definitions fl3.43p flA.Sp . we readily infer the 
estimates 


40) 


L(no) 




(A.37) 


?(i) ^ 4 m+i(^'o) 


l^*,M+i(^o,r)| ^ ^M+l)\ [4(^+^+2,nor )lnr -0(M+A+2,nor ) 
these, along with the bounds in Eq.s flA.lip (IA.20p . imply in turn 


2 \ 2 (Af+A+ 2 ) . 


fil%+i(vo,r)\ « (^^qXFT1^4t+i(»'o)l liir" + 


-2(Ar+A+2) 


r 


2 , (^+1)! 


(A.38) 


(M+A+2)^ 


2(AfH-AH-2) 


14 


For example, we have the following chain of equalities: 


pvor^ 

Inz; = r- 2 (’"+^+i) / dw \niw/r^) = 

Jo 


= r“2(™+“+i)/ dwe-^w'^+^ 


In w—In R 


^ ^- 2 (m+A+l) 


g(TO+A+l, Vo r^)— 7 (m+A+l, nor^) In R 
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Summing up, Eq.s fl3.4ip fl3.42|) follow from Eq. flA.291) and the relations in Eq.s 
(1A.34D (1A.35P with M = K + i — I; the remainder estimate fl3.44p descend easily 
from Eq.s flA.31[) (lA.38p . which also give explicit expressions for the constants in 
fl3.44p . Finally, using the asymptotic expansions in Eq.s flA.lip flA.23p (respectively 
holding for the lower incomplete gamma and for the function g( , )), Eq.s flA.34p 
flA.3811 imply Eq.s fl3.45p fl3.46p . 


B Appendix. An alternative representation for 
the bulk energy 

In subsection l3.8l we obtained an integral representation for the renormalized bulk en¬ 
ergy (see Eq. 03.551) ) to be evaluated numerically (Q). In the present appendix 
we derive an alternative representation for the regularized energy allowing to 
express it in terms of the Riemann zeta function (see [20], p. 602, Eq. 25.2.1) 

- 1-00 .. 

C(s) := — for s G C with > 1 . (B.l) 

n=l 

Next, the renormalized energy is computed via the zeta approach, using the 
well-known analytic continuation of ^ to the whole complex plane. 

Let us stress that the methods discussed here only works if the background potential 
is isotropic. Nonetheless, they can be of some interest since, for example, they allow 
us to perform a direct comparison with the results of Actor and Bender |2] in the 
case with d = 3 (see subsection IB.3|1 . 

In the forthcoming subsection IB. ll we introduce a family of integrals and discuss some 
relations allowing to express them in terms of the Riemann zeta; in the following 
subsection IB. 2 1 the regularized bulk energy is represented in terms of these integrals 
and its analytic continuation is discussed. In the concluding subsection IB. 31 we use 
the results of subsections IB. 1 1 [B. 2 1 to derive, as examples, the renormalized energy 
Eren cases with d G {1, 2, 3} . 

B.l The integrals Jni^) a recursive relation. Let us consider the family 
of integrals 

^- 1-00 

J„(s) := / dr —for nG {1, 2, 3,...}, sG C with 3fJs > u (B.2) 
Jo sinh r 

(of course, the restriction on 3fJs arises in order to guarantee the convergence of the 
integral J„). 

pointed out in the cited subsection that the main advantage of this approach is that it 
work also for configurations where the background harmonic potential is not isotropic. 


42 














































First note that, for n = 1 and n = 2, the above integrals are known to be related to 
the Riemann zeta function more precisely, we have (see [20], p.604, Eq.s 25.5.8 
and 25.5.9) 

ai(s) = 2(1-2-*) r(s)C(s) , (B.3) 

a2(s) = 2-(*-2)r(s)C(s-l) . (B.4) 

Next, let us £x n G {1, 2, 3,...} and consider the elementary identity 


(f ( 1 


- 


n{n+l) n 


+ 


dr'^ V sinh"'r J sinh’^’'"^r sinh'^r ’ 
using this result and the definition flB.2p . we infer (for s G C with 3fJs > n + 2) 


(B.6) 


fln+2(s) — — 


n 


n + 1 




n(n+l) Jq 


+00 j2 

dTT^-^ ^ 


V sinh"'r 


(B.6) 


Concerning the integral appearing in the right-hand side of the above equation, 
integrating by parts two times we obtain 




'0 


V sinh"'r 


'‘+CO 


(s-l)(s-2) 


nS— 3 


(B,7) 


dr 


sinh"r 


= (s-l)(s-2)a„(s-2) . 


Summing up, Eq.s fIB.Op flB.7p give the recursive relation 

a„+2(s) = 


^ J„(s) + ^ a„(s - 2) (s G C, 3fis > n + 2) . (B.8) 


n-|-l 


n(n-M) 


This result and the explicit expressions flB.3p (lB.4p allow to express any integral 
Dfn(s) in terms of (suitable linear combinations of) the Riemann zeta C,] needless to 
say, the analytic continuation of Jn is then determined by the well-known analytic 
continuation of C,. 

B.2 The bulk energy in terms of the integrals Using Eq. fl3.53p for 

the regularized energy 77“ with n = 0 and the explicit representation fl3.52p for the 
heat trace K{t), we obtain 


= 


K 


^+oo 


dt 


n —3 

t 2 


2d+iY{^) sinh'^(fc2t) ’ 


(B.9) 


making the change of variable r := k^i G (0, -|-cx)) and recalling the definition flB.2p . 
we infer 


= 


k 


2^+1 r(H^) \k) 


T 2d 


u 


(B.IO) 
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The recursive relation fIB.SD for 3^ and the identities fIB.Sp flB.dp allow us to deter¬ 
mine the analytic continuation of to the whole complex plane via Eq. fIB.lOp . 
The renormalized bulk energy is defined according to Eq. fl2.39p . setting 


Even Ep 




u=0 


or := E" 


u=0 


The hnal expression of E^^'^ for arbitrary spatial dimension d is too complicate to 
be reported here; in the next subsection we consider, as examples, the cases with 
d e {1,2,3}. 

B.3 The bulk energy in terms of the Riemann zeta function for d G {1, 2, 3} 

As mentioned at the end of the previous subsection, here we exemplify the general 
methods developed in the previous subsections IB. lIlB.21 for the computation of 
in the cases with d G (1, 2, 3}; we also show that the results found are in agreement 
with those derived in Section 0] using another approach. 

The case d = 1. Using Eq.s flB.3p fIB.lOp and setting u = 0, we infer 


Even ^ 




(B.ll) 


Evaluating numerically this expression, we have (in agreement with Eq. fl4.15lB 


E^^^ = fc(0.0430546468...) . 

The case d = 2. Eq.s (lB.4p and (IB.lOp with u = 0 imply 

v/2 A 2 

Numerical evaluation of the above result gives 

Even ^ _fc(o.0180207590...) , 


(B.12) 


(B.13) 


(B.14) 


which agrees with the one of Eq. fl4.33l) (the last digit being different due to trun¬ 
cation approximation). 

The case d = 3. This time we have to resort to the recursive relation flB.8p (here 
employed with u = 1) as well as to the identity flB.3l) . Then, using once more Eq. 
fIB.lOl) with M = 0, we infer 


Even ^ ^ 




(B.15) 


We state that the above result coincides with the one of Eq. (4.4) in [2], which in 
our language reads (see footnote [12]) 


^ren 




1 



(B.16) 
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where C( 5 ) indicates the analytic continuation of the Hurwitz zeta function; this 
statement follows is easily proved using the known identities (see [20], p. 607, Eq. 
25.11.3 and p. 608, Eq. 25.11.11) 


C(s,a+ 1) = C(s, a) - a 


C(<.,-|=(2‘-1)CW. 


(B.17) 

(B.18) 


Finally, let us note that evaluating numerically the expression in the right-hand side 
of Eq. (1B.15P we have 

Even ^ _fc(o.0078607118...) , (B.19) 

in agreement with Eq. (I4.52p . 
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